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The  characteristic  functional  of  the  derivative  ^(t)  of  a  Markov  process  <ti(t)  and 
the  related  multiplicative  process  cr(t ),  which  obeys  the  stochastic  differential  equation 
io(t)  =  (A  +  <i>(t)B)o (t) ,  have  been  studied.  Exact  equations  for  the  marginal  character¬ 
istic  functional  and  the  marginal  average  of  o(t)  are  derived.  The  first  equation  is 
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background.  The  equation  of  motion  for  the  marginal  average  of  a(t)  is  solved  for  various 
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I.  INTRODUCTION 


The  equation  of  motion  for  the  density  operator  of  an  atom  in  a  finite- 

bandwidth  laser  field  or  the  equation  for  the  regression  of  the  atomic 
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dipole  correlations  assumes  the  general  form  * 

i  ^  =  (A  +  ♦(t)B)o  ,  (1.1) 

where  A  and  B  are  linear  operators  in  Liouville  space,  which  act  on  the 
Liouville  vector  o(t).  Here  ♦(t)  represents  the  laser  phase,  which  is 
considered  to  be  a  real-valued  stochastic  process.  The  fluctutating  phase 
broalens  the  laser  line,  but  the  atom  responds  to  the  instantaneous 
frequency  shift  ♦(t),  which  is  the  time-derivative  of  the  laser  phase.  The 
process  (^(t)  is  again  a  stochastic  process,  and  via  Eq.  (1.1)  the  state  of 
the  atom  or  the  correlation  functions  o(t)  become  stochastic  quantities. 

The  issue  in  quantum  optics  is  then  to  solve  the  multiplicative  stochastic 
differential  equation  (1.1)  for  the  average  <o(t)>.  The  first  solution  was 
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obtained  by  Fox,  who  assumed  the  process  ^(t)  to  be  Gaussian  white  noise, 

which  corresponds  to  a  diffusive  Gaussian  phase  ^(t)  (the  Wiener-LAvy 

process).  This  result  was  generalized  to  a  Gaussian  process  A(t)  with  a 

finite  correlation  time  and  an  exponentially  decaying  correlation 

function^  ^  (the  Ornstein-Uhlenbeck  process),  and  to  a  process  ACt),  which 
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is  again  diffusive,  but  not  Gaussian  ’  (the  independent- increment  process). 
Furthermore,  Eq.  (1.1)  can  be  solved  for  <o(t)>  if  we  have  A(t)  as  a  Markov 
random-jump  process, which  models  a  multimode  laser. 

In  these  examples  the  solvability  of  the  problem  relies  on  the  Gaussian 


property  of  ♦(t),  or  hinges  on  the  prescribed  stochastics  of  A(t).  This 
implies  that  the  process  A(t)  is  actually  considered  to  be  the  driving 


process.  For  a  single-mode  laser  In  general,  however,  the  phase 
fluctuations  ^(t)  are  specified  rather  than  the  derivative  ^(t)  of  this 
process.  A  prime  example  would  be  the  atomic  response  to  phase- locked 
radiation, as  it  is  generated  for  instance  by  some  ring  lasers. In  this 
paper  we  shall  develop  a  general  method  to  solve  Eq.  (1.1)  for  the  case  that 
^(t)  is  a  given  Markov  process.  The  formal  theory  will  be  exemplified  by  a 
specific  choice  for  ^(t),  which  models  phase-locked  radiation.  Furthermore, 
we  shall  study  the  time-derivative  of  ♦(t)  itself  and  extract  the 
stochastics  of  ^(t)  from  the  properties  of  ^(t). 


II.  THE  STOCHASTICS  OF  (»(t) 

18 

Let  us  define  the  phase  ♦(t)  as  a  homogeneous  Markov  process.  Then 

its  stoch  sties  is  fixed  by  the  probability  distribution  P(^,t)  and  the 

conditional  probability  distribution  (“tiO),  which  has  the 

significance  of  the  probability  density  for  the  occurrence  of  ♦(t+x)  *  (>2  if 

^(t)  >  For  a  homogeneous  process  this  is  independent  of  t  by 

definition.  The  higher-order  statistics  is  now  determined  by  the  Markov 
19 

property.  From  the  obvious  relation 


p.  .  (♦lo' )?(♦', t  ) 


p(^,t). 


t  >  t. 


(2.1) 


it  follows  that  it  is  sufficient  to  prescribe  the  probability  distribution 
P(^,t)  for  a  single  time-point  t^  only.  The  time-evolution  towards  t  >  t^ 
can  then  be  found  from  Eq.  (2.1)  and  P  (^|^'). 

0  jg 

The  conditional  probability  distribution  obeys  the  Master  equation 


—  P^(^3|^p  -  Jd^2  {W(^3|^2)  - 


(2.2) 


with  W(^'l^)  >  0  as  the  transition  rate  of  the  process  from  ^  to  and 


4 


a(^)  -  /<!♦'  W(^'  1^)  , 


(2.3) 


which  is  the  loss-rate  of  independent  of  the  final  value  .  The  Initial 


condition  for  Eq.  (2.2)  reads 
Po(^3ki)  ■  ^(♦3-^i)  . 


(2.4) 


so  a  given  W(^' |^)  determines  P^(^3l^2^  every  x  >  0.  Hence  the 
stochastics  of  a  homogeneous  Markov  process  <^(t)  is  fixed  as  soon  as  P(^,tQ) 
and  W(^'|4i)  are  prescribed.  These  functions  will  from  now  on  be  assumed  to 
be  given. 


III.  THE  CHARACTERISTIC  FUNCTIONAL 


A  convenient  way  to  represent  the  stochastic  properties  of  a  stochastic 

9  20 

process  is  by  means  of  its  characteristic  functional.  '  Since  we  are 
concerned  with  the  process  ♦(t),  we  define 


Z^lk]  -  <exp(-i/ds  ♦(s)k(s))>  , 


t  >  t. 


(3.1) 


which  is  a  functional  of  the  test  function  k(t).  Here  the  angle  brackets 

denote  an  average  over  the  stochastic  process  ♦(t)  or  ♦(t),  whatever  is 

prescribed.  A  general  method  to  evaluate  Z^lk]  for  the  case  where  ♦(t)  is  a 

21 

homogeneous  Markov  process  has  been  given  by  van  Kampen. 

Knowledge  of  the  characteristic  functional  Z^[k)  determines  completely 
the  stochastics  of  ^(t),  which  can  be  seen  as  follows.  Choose  k(s)  as  the 
sequence  of  6- functions 


k(s)  ”  ' 


^  >  ^0  * 


(3.2) 


and  take  t  ■  *  in  (3.1).  Then  we  f ind 


/  .V.*  - 


Z^[k]  -  <exp(ik  iCt  )  +  ...  +  ik  ♦(!,))>  ,  (3.3) 

•  n  n  11 

which  is  the  moment -general ins  function  of  ^(1).  If  we  write  z  (k  ,t  ; 

n  n  n 

...  ;k,  ,t,}  for  Z  [k],  then  we  can  obtain  the  moments  of  it»(t)  =  A(t) 

11  • 

according  to 


<»(t  )  ...  ♦(tj)>  -  (-1)"  ...  !  ...ik  .t  >1 

n  1  k»...»k,  *0 

n  1 


and  the  probability  distributions  by 


(3.4) 


n  n  n  i  l 


(2ii) 


.dk  ...  dk, 
n  '  n  1 


-ik  ♦  -  ...  -ik,*, 

*  e  z  (k  ,t  ;  ...  ;k,,t,)  , 

n  n  n  11 


(3.5) 


where  we  have  introduced  in  order  to  distinguish  from  the  probability 


distributions  for  *(t)  itself. 


IV.  THE  MARGINAL  AVERAGE 
A.  GENERAL 

The  exponential  in  Eq.  (3.1)  is  a  functional  of  both  k(t)  and  *(t),  so 

it  depends  on  the  values  of  *(t)  in  the  complete  interval  [t^.t).  After  the 

average  has  been  taken  it  will  be  only  a  functional  of  k(t).  The  general 

attempt  to  evaluate  averages  of  a  functional  is  to  derive  an  equation  for 

the  average.  For  subsequently  solving  this  equation  for  functionals  which 

involve  Markov  processes,  this  scheme  is  most  conveniently  carried  out  by  an 
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intermediate  Introduction  of  Burshtein's  marginal  averages.  Since  in  our 
problem  the  stochastics  of  *(t)  is  assumed  to  be  given,  the  appropriate 
marginal  characteristic  functional,  which  is  related  to  Z^[k],  should  be 


defined  as 


(4.1) 


-  <6(^(t)-#Q)  exp(-i  Jds  ♦(s)k(s))>  , 

for  t  >  t^.  The  initial  value  is  then 

Qt^[^0»*^l  -  <6(^(tQ)-^Q)>  -  P(^o‘*^o)  *  (^-2) 

and  Z^[k]  follows  from  Q^[^Q,k]  according  to 

Z^.[k]  *  Jd^Q  Qj.[«Q.kl  .  (4.3) 

For  t  *  t^  we  find  with  Eq.  (4.2)  [k]  =  Jd^Q  P(^Q,tQ)  =  1 ,  in  agreement 

with  Eq.  (3.1). 

In  order  to  derive  an  equation  for  the  time-evolution  of  Q^[^Q,k], 
we  first  increase  t  by  a  small  amount  At  >  0.  This  gives 


<6((>(t+At)-(>Q)  exp{-i((>(t+At)-(>(t))k(t)} 
t 

*  exp(-i  Jds  ^(s)k(s))> 


(4.4) 


Subsequently,  we  expand  the  exponential  functional  of  ^(s)  in  a  series,  and 
we  take  the  average  in  (4.4)  term  by  term.  Hereafter,  we  apply  the  Master 
equation  (2.2)  for  P^+At^^^^O^  take  the  limit  At  0.  This  yields  an 

equation  for  the  marginal  average,  and  explicitly  we  find 

at  (^(♦ol^^  ■  a((>)6((>Q-(>)} 

-i(i>--^)k(t) 

X  e  Qj,U,k3  .  (4.5) 


The  Markov  process  (>(t)  is  characterized  by  P(^Q,tQ)  and  W(^q|^),  which 
respectively  determine  the  initial  value  and  the  time-evolution  of  Q^[^q 


.k) 


For  a  specific  choice  of  W(^q|^),  we  have  to  solve  Eq.  (4.5),  after  which 
the  characteristic  functional  Z^[k]  can  be  obtained  from  Eq.  (4.3). 

Notice  the  resemblance  between  the  result  (4.5)  and  the  Master  equation 
(2.2).  If  we  multiply  Eq.  (2.2)  by  P(^j,tQ),  take  x  ■  t  -  t^  and  apply  the 
relation  (2.1),  we  find 

^  P(^Q,t)  =  Jd^  {W(*q|^)  -  a(^)6(^Q-^)}  P(#,t)  ,  (4.6) 

which  is  the  Master  equation  for  P(^Q,t).  This  equation  is  identical  to  Eq. 
(4.5),  including  the  initial  condition  (4.2),  if  we  set  k(t)  =0.  On  the 
other  hand,  it  follows  from  Eq.  (4.1)  that  Q^[^Q,k]  =  <6(^(t)-^Q)>  =  P(^Q,t) 
if  we  take  k(t)  =  0  ,  so  that  in  this  case  Eq.  (4.5)  should  indeed  reduce  to 
Eq .  (4.6). 

B.  INDEPENDENT  INCREMENTS 

In  order  to  display  the  usefulness  and  applicability  of  the  marginal- 
functional  approach,  we  consider  an  example.  Let  us  specify  the  transition 
rate  by 

W((>q|<>)  =  T  w((>Q-^)  ,  y  >  0  ,  (4.7) 

where  the  function  w(n)  is  normalized  as 

Jdn  w(n)  *  1  .  (4.8) 

The  stochastic  process  ♦(t)  will  be  defined  on  the  real  axis,  with  -“  <  ^  < 
•.  The  assertion  (4.7)  states  that  the  probability  for  a  transition  (>  -♦ 
depends  only  on  the  phase  difference  ♦q  -  and  from  Eq.  (2.3)  we  find  that 
a(^)  ■  Y,  so  that  the  total  loss-rate  for  ^  is  Independent  of  This  is  a 
diffusion  process,  and  it  is  commonly  referred  to  as  the  independent- 
increment  process.  As  an  Initial  condition  for  the  probability 
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distribution,  we  take 

P(0,tQ)  =  6(^)  .  (A. 9) 

Comparison  of  the  Master  equations  for  and  P(^,t)  then  shows  that 

the  probability  distribution  and  the  conditional  probability  distribution 
are  related  according  to 


(4.10) 


The  Master  equation  (4.6)  for  P(^,t)  can  be  solved  by  Fourier 
transformation  with  respect  to  If  we  write 


f(p 


,t)  =  f  M  . 


(4.11) 


which  has  f(p,tQ)  ®  1  as  the  initial  condition,  then  the  solution  of  Eq. 
(4.6)  is  immediately  seen  to  be 

Y(w(p)-l)(t-t_) 

P(p,t)  =  e  ,  t  >  tp  ,  (4.12) 

in  terms  of  the  Fourier  transform  w(p)  of  w(^).  Note  that  w(0)  =  1,  as  a 
result  of  the  normalization  (4.7).  Along  the  very  same  lines  we  can  solve 
Eq.  (4.5)  for  the  Fourier  transform  0^(p»k].  We  obtain 

t  « 

Oj^Ip.k]  »  exp(-Y  /  ds  J  d^  (1-e^*^^^  ^ )w((>) ) 

t 

0 

after  which  the  characteristic  functional  follows  from 
Z^.[k]  -  0^(0, k]  , 

g 

which  yields  the  familiar  result  . 


(4.14) 
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V.  THE  MARGINAL  MOMENTS 
A.  GENERAL 

If  we  take  k(s)  as  the  sequence  of  delta  functions  (3.2)  in  the 
definition  (4.1)  of  the  marginal  characteristic  functional,  it  assumes  the 
form 


n 

Q^.l^Q.k]  =  <6((|i(t)-(j)Q)exp(i  ^  kjj^  0(t-tjj^))>  ,  (5.1) 

i=l 

with  ♦(t)  =  ♦(t)  and  9(t)  the  unit-step  function.  Just  as  we  can  find  the 

moments  <4)(t  )  ...  ii(t,)>  of  <Ji(t)  from  Z  [k],  we  can  obtain  the  marginal 
n  1  * 

moments  <6((>(t)-^Q)4)(t^)  ...  ip(tj)>  from  Qj.[^Q,k].  Obviously  the  integral 
over  (>Q  of  the  marginal  moments  yields  the  moments.  The  characteristic 
functional  Z^[k]  becomes  independent  of  t  if  t  >  t^^  for  all  I,  but  Q^[i>Q,k] 
remains  t ime- dependent .  This  is  due  to  the  appearance  of  6((>(t)-(>Q) . 
Furthermore,  the  time  t  is  a  dynamical  variable  in  Eq.  (4.5),  so  that  care 
should  be  exercised  in  the  time-ordering.  The  marginal  moments  follow  from 
by  differentiation,  according  to 


<6(^(t)-(>Q)4»(t^)  ...  <(>(t^)>  0(t-t^)  ...  0(t-t^)  * 

3k  3k  = 

n  1  n 


.  =  k^  =  0 


.  (5.2) 


Equation  (4.5)  for  Q^('>Q,k]  implies  an  equation  for  the  marginal 
moments.  First,  we  note  that 

exp{-i((>Q-<>)k(t)}Q^[(>,k]  = 

n 

<6((>(t)-(>)exp(  i  ^  k^{(<|)Q-(>)6(t-t^)  +  '<'(t^)0(t-tj^)})>  . 

i=l 


(5.3) 
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After  substituting  this  expression  in  the  right-hand  side  of  Eq.  (4.5), 
differentiating  with  respect  to  k^,  ...  .k^,  setting  ®  •••  ~  ® 

integrating  over  time,  we  obtain 

<6(4»(t)-(>Q)i(i(t^)  ...  e(t-t^)  ...  9(t-tj)  =  Jd(>  {W((>q|()i)  -  a((|))6((tiQ-<>)} 

t 

X  J  dt'  <6(i>(t' )-<>){(^fv-'>)5(t' -t  )  +  4i(t  )e(t'-t  )} 

'  U  n  n  n 

0 

...  {(((.Q-().)6(t'-tj)  +  t|.(t^)0(t'-t^)}>  .  (5.4) 

Vfhen  we  set  t  >  t^^  for  all  I,  we  have  a  Master-like  equation  for  <5((>(t)-(>Q) 

it)(t  )  ...  il)(t,)>,  and  the  lower-order  marginal  moments  <6((>(t)-()i- )4)(t  )  ... 
n  i  Urn 

'>'(t^)>  with  m  <  n  appear  as  inhoroogenous  terms.  Hence,  Eq.  (5.4)  should  be 

solved  successively  for  n  =  1,  n  =  2 .  We  note  that  Eq.  (5.4) 

provides  an  explicit  expression  for  <i(»(t^)  ...  ij)(t^)>  in  terms  of  the  lower- 
order  marginal  moments  after  an  integration  over  4>q.  Indeed,  from  the 
property 

/d(>Q  {W(^qI^)  -  a(0)6(^Q-0)}  =  0  ,  (5.5) 

the  term  with  <6(i>(t' )-0)4/(t^)  ...  il)(tj^)>  on  the  right-hand  side  of  Eq. 

(5.4)  vanishes  after  an  integration  over  (>  . 

u 

B.  LOWEST  ORDERS 

In  order  to  exhibit  clearly  the  structure  of  the  equation  for  the 
marginal  moments,  we  consider  the  cases  n  =  1  and  n  =  2  in  some  more  detail. 
After  a  slight  rearrangement,  Eq.  (5.4)  for  n  =  1  can  be  written  as 

<5(0(t)-*Q)ii>(t^)>  =  fd<p  {Vl((p^l(p)  -  a((>)6(4.Q-<i)} 

t 

*  {(4in-(t)P(<|i,t  )  +  Jdt'  <6((|>(t' )-<(i)ii'(t,  )>}  , 

^  ^  tj  ^  (5.6) 


for  t  2  t^.  This  integral  equation  in  time  is  equivalent  to  the 
differential  equation 


<6(«(t)-4Q)*(t^)>  -  Jd^  {W(^qI^)  -  a(^)6(^Q-^)} 

X  <6(^(t)-^H(tj)>  ,  (5.7) 


together  with  the  initial  condition 

<6((>(t^)-^Q)ip(t^)>  =  Jd^  {W(^qI^)  -  a((>)6((jQ-(»)}(^Q-(>)P(*,t^)  .  (5.8) 

The  equation  for  the  first  marginal  average  <6(^(t)-^Q)4i(tj)>  is  identical 
to  the  Master  equation  (2.2),  but  with  a  different  initial  value. 

Integration  of  (5.8)  over  ♦q  yields 

«i(tp>  =  Jd(>Jd^Q  W(^Q|^)(^Q-^)P(^,tj)  ,  (5.9) 

which  expresses  explicitly  the  average  of  <i(»(tj^)>  in  the  given  functions 
W(^qI^)  and  P(^,tj).  With  the  aid  of  the  Master  equation,  we  can  cast  (5.9) 
in  the  form 


<*(t,)>  =  Jd<>  (K  ^  P(«,t,) 


3t 


1 


(5.10) 


as  it  should  be. 

The  solution  of  Eq.  (5.6)  for  <6(^(t)-^Q)^(tj)>  provides  the  input  for 
the  explicit  expression  for  the  two-time  correlation  function,  which  becomes 

<4«(t2)i(»(t^)>  =•  |d(>/d^Q  {W(<>q|^)  -  a(^)6(^Q-(>)} 

X  {((>Q-^)^6(t2-t^)P(<>,t2) 

+  (♦Q-(>)(<6(t(t2)-^)*(tj)>e(t2-tp 

+  <6(^(t,)-0)(|)(t_)>e(t,-t-))} 


(5.11) 


The  appearance  of  shows  that  the  time-derivative  of  any  Markov 

process  is  6-correlated  with  a  continuous  background. 

C.  RANDOM  JUMPS 

Equation  (5.11)  for  instance  might  seem  awkward,  but  it  is  really 
straightforward  in  its  application.  Let  us  illustrate  this  with  an  example. 
Consider  the  random- jump  process  <^(t),  defined  as  a  stationary  process  with 
transition  rate 

«(♦!♦')  =  Y  ?(♦)  .  Y  >  0  ,  (5.12) 

in  terms  of  an  arbitrary  probability  distribution  P(^).  Eq.  (5.12)  is 

equivalent  to  the  statement  that  the  probability  for  a  transition  -►  (>  is 

13 

Independent  of  the  initial  value  .  From  Eq.  (5.9)  we  immediately  derive 


<*(tj)>  =  0  ,  (5.13) 

which  is,  in  view  of  (5.10),  necessary  for  a  stationary  process.  From  (2.3) 
we  obtain  a((>)  »  y,  and  the  solution  of  Eq.  (5.7),  with  initial  value  (5.8), 
is  readily  found  to  be 


-Y(t-t, ) 

<6(^(t)-(>g)i|p(t^)>  «  yP((>Q)((>Q-bj)  e 


Here  we  have  introduced  the  moments  of  P(^)  as 


t  >  t^  .  (5.14) 


b^  *  Jd^  P((>)  , 

which  are  parameters  of  the  process  <>(t).  Solution  (5.14)  can  be 
substituted  into  Eq.  (5.11),  which  gives  the  correlation  function 

2  -Y)t--t  I 

<*(tj)*(t2)>  =•  Y(b2-bp{26(t^-t2)  -ye  M  . 

for  all  tp  t2. 


(5.15) 


(5.16) 


From  (5.15)  it  follows  that 
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b^-b^  >  0  ,  (5.17) 

so  that  for  t^  *  t^  the  correlation  (5.16)  is  negative.  For  t^^  =  t^  the  4- 

2 

function  dominates  the  negative  term,  so  that  <))>(tj^)  >  is  positive,  as  it 
should  be. 


VI.  THE  MULTIPLICATIVE  PROCESS 

So  far  we  have  considered  the  stochastics  of  ^(t)  itself,  and  its 
characteristic  fxmctional.  In  this  section  we  shall  generalize  the  method, 
in  order  to  solve  the  multiplicative  equation  (1.1).  To  this  end  we  write 
the  formal  solution  of  (1.1)  for  the  stochastic  vector  o(t)  as 


-iA(t-t-)  C 

ar(t)  =  e  T  exp[-ijds  (>(s)8(s) ]o(t-)  , 

^0 

where  T  is  the  time-ordering  operator  and  B(t)  is  defined  as 


(6.1) 


iA(t-t^)  -iA(t-t.) 
B(t)  =  e  °  B  e  ^ 


(6.2) 


In  close  analogy  to  the  definition  of  f*'  E*!-  (^*1),  we  now 

introduce  the  marginal  average  of  o(t)  by 


;(^Q,t)  -  <6((>(t)-«Q)o(t)>  .  (6.3) 

Then  we  substitute  the  expression  (6,1)  for  o(t)  and  replace  t  by  t  +  At, 
which  gives  a  similar  formula  as  Eq.  (4.4).  That  this  can  also  be  done  for 
the  time-ordered  exponential  is  sometimes  referred  to  as  the  semi-group 
property  of  the  evolution  operator.  Along  the  same  lines  that  led  to  Eq. 
(4.5)  we  now  f ind 


^3t  *  AC(^Q,t) 


f  "i(*>Q“^)B 

+  ijdi>  {W((>  |(>)  -  a(^)<5((>--t)}e  ^(♦,t) 


(6.4) 
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or  equivalently 

3  f 

(i^  -  A  +  ia((>Q))  C((>Q,t)  =  ijd<>  W((>Q|(>)e  ^(ip,t)  .  (6.5) 

Notice  that  the  operator  B  appears  in  the  exponential,  rather  than  S(t), 
as  could  be  expected  by  analogy  with  the  characteristic  functional.  For  a 
given  stochastic  process  ^(t),  e.g.,  a  given  and  PC^.t^),  we  have  to 

solve  Eq.  (6.4)  with  the  initial  condition 

C(^o,to)  =  <6(^(tQ)-^Q)o(tQ)>  ,  (6.6) 

after  which  <o(t)>  follows  from 

<o(t)>  =  |d.>Q  C((>Q.t)  .  (6.7) 

For  a  given  non-stochastic  state  0(tQ),  the  initial  condition  reduces  to 

C(«>Q,tQ)  *  P(^Q,tQ)  o(tQ)  ,  (6.8) 

which  differs  from  (6.6)  by  the  fact  that  there  are  no  initial  correlations. 
This  means  that  the  process  a(t)  has  no  memory  to  times  smaller  than  t^,  and 
consequently  its  evolution  for  t  >  t^  is  completely  determined  by  its 
initial  state  o(tQ).  It  was  emphasized  by  Arnoldus  and  Nienhuls  that  the 
comnon  choice  ^(♦g»tg)  =  P(^Q,tQ)  <o(tp)>  is  merely  an  approximation  which 

only  holds  for  small  correlation  times  of  ♦(!). 

VII.  SOLUTIONS 
A.  INDEPENDENT  INCREMENTS 

Equation  (6.5)  for  the  marginal  average  of  o(t)  can  be  solved  for  the 
independent- increment  process  with  the  same  procedure  as  in  Section  IV, 
where  we  obtained  the  characteristic  functional.  If  we  adopt  the  Fourier 


transform 
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C(p.t)  =  C(*.t)  =■  ,  (7.1) 

•00 

where  the  second  equality  follows  after  application  of  Eq.  (6.3),  then 
<a(t)>  can  be  found  from 

<o(t)>  -  C(0,t)  .  (7.2) 

With  the  technique  of  Section  IV  we  can  find  C(p,t),  and  if  we  differentiate 
the  result  with  respect  to  time,  we  find 

i^  C(p,t)  =  (A  -  iW(p))C(p,t)  ,  (7.3) 


with 

ft(p)  =  Y/**dn  ®^)  w(n)  .  (7.4) 

•00 

The  operator  ft(p)  accounts  for  the  phase  fluctuations.  If  we  set  p  =  0  in 
Eq.  (7.3),  we  achieve  the  equation  for  <o(t)>,  with  solution 

-i(A-ift(0))(t-t-) 

<o(t)>  =  e  ^  <o(tQ)>  ,  (7.5) 

for  t  >  tQ.  We  note  that  <o(t)>  can  be  expressed  in  terms  of  <o(tQ)>  for 
this  process,  so  that  there  are  no  initial  correlations  for  the  diffusion 
process.  The  process  ♦(t)  has  no  memory,  and  with  the  results  of  Section  V 
it  can  be  shown  that  ♦(t)  is  indeed  delta-correlated.  This  means  that 
<<>(t^)  ...  (>(tj^)>  for  all  n  contains  only  6-functions,  which  implies  the 
factorization  in  (7.5). 

A  special  case  arises  if  we  take 


yw(n)  ■  Y  ^(n)  +  X  6"(n)  ,  x  >  0 


(7.6) 
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where  the  primes  on  the  6-function  denote  differentiation  with  respect  to 
its  argument.  It  is  easy  to  check  that  this  process  is  the  Wiener-L6vy 
process,  or  the  phase-diffusion  process.  If  we  substitute  (7.6)  into 
(4.12),  we  find  that  P(^,t)  is  Gaussian,  and  obviously  this  is  the  only 
Gaussian  limit  of  the  diffusion  process.  The  operator  0(p)  in  (7.4)  reduces 


ft(p)  =  X(p-B)  , 


and  the  equation  for  <a(t)>  becomes 


i-j7  <o(t)>  =  (A-iAB^)<a(t)>  , 

at 


(7.7) 


(7.8) 


which  is  the  result  of  Fox  . 

B.  QRNSTEIN-UHLENBECK  PROCESS 

The  diffusion  process  has  no  memory  and  is  essentially  non-stationary. 
The  initial  distribution  P(^,tQ)  »  6(^)  diffuses  over  the  whole  (^-axis, 
-«<♦<<».  The  inclusion  of  a  finite  memory-time  can  stabilize  this 
process.  Let  us  define  the  transition  probability  as 


W(^qI^)  -  a(^)6(^Q-0)  »  A  6"(<>q-^)  +  y  ♦  6'('>o-^)  . 


A  >  0  , 

Y  >  0  . 


(7.9) 


Then  the  Master  equation  (4.6)  for  P(i>,t)  becomes  the  Fokker- Planck 
equation 


^  p((>,t)  =  (A  af  P(<>.t)  , 

3^ 


which  has  the  solution  for  t  -►  «> 


(7.10) 


P(^)  -  (2iro^)'* 


o  »  A/y 


(7.11) 


•  ^  k  _ii  _  ^ 


This  ?(♦),  together  with  W(^q|^)  from  (7.9),  defines  a  stationary  Gaussian 

Markov  process,  the  Ornstein-Uhlenbeck  process.  In  the  limit  -y  -+  0  and  X 

finite  (so  o  ->•“),  the  process  ♦(t)  reduces  to  the  Wiener-L6vy  process  from 

the  previous  paragraph.  From  (7.11)  we  see  that  ^(t)  is  centered  around  ^  » 

2 

0.  The  distribution  is  Gaussian  with  a  variance  a  around  the  average  ^ 

0.  The  preference  for  ^  >  0  expresses  that  this  process  can  be  considered 
as  a  model  for  phase- locked  radiation. 

With  the  specific  choice  (7.9)  for  the  transition  rate,  the  Master 
equation  (6.5)  assumes  the  form  of  a  second-order  partial  differential 
equation.  We  obtain 

(ig|  -  A  +  iAB^)  ^(♦.t)  *  iY{o^(2iB  +  ’ 

(7.12) 

2 

In  the  limit  y  0  and  A  ■  yo  finite,  we  recover  (the  Fourier  inverse  of) 
Eq.  (7.3)  with  ft(p)  from  Eq.  (7.7). 

In  order  to  obtain  a  solution  of  Eq.  (7.12),  we  start  with  a  Fourier 
transform  with  respect  to  The  transformed  equation  then  reads 

(i^  -  A  +  iAB^)C(p,t)  “  -iY{o^(p-2B)p  +  (p-B)  ^}C(p,t)  ,  (7.13) 


which  is  still  a  partial  differential  equation.  Since  we  are  interested  in 
C(0,t)  *  <o(t)>,  the  obvious  approach  would  be  a  Taylor  expansion  around 
p  *  0.  This  yields  however  a  cumbersome  inhomogeneous  four-term  recurrence 
relation  for  the  Taylor  coefficients.  This  can  be  avoided  by  the 
transformation^^ 


S(p,t)  -  P(p)  g(p,t) 


(7.14) 


which  defines  g(p,t).  The  Fourier  transform  of  the  probability  distribution 
is  explicitly 


P(p)  .  P  .  (7.15) 

and  in  particular  we  have  f(0)  =  1.  The  equation  for  g(p,t)  becomes 

(igl-  -  A  +  iAB^)g(p,t)  »  -iY{p  *  B(o^p  +  ^)}g(p,t)  ,  (7.16) 

and  it  has  to  be  solved  for 

g(0,t)  “  <a(t)>  .  (7.17) 

Let  us  define  the  Taylor  coefficients  by  the  expansion 

00 

g(p,t)  =  *  (7.18) 

n=0 

which  can  be  inverted  as 

X  (t)  «  <o(t)  (^)"  •  (7.19) 

n  oip  p=0 

Substitution  of  (7.18)  into  (7.16)  then  gives  the  equation  for  the  Taylor 
coefficients 

(i^  -  A  +  iAB^  +  iYn)x  (t)  =  YB(no^  x  (t)  -  x  .,(t))  ,  (7.20). 

ot  n  n-l  n+l 

which  has  to  be  solved  for 


XQ(t)  »  <o(t)>  .  (7.21) 

Equation  (7.20)  looks  like  a  homogeneous  three- term  recurrence  relation,  but 
it  will  be  shown  below  that  the  time -derivative  d/9t  gives  rise  to  an 
inhomogeneous  contribution.  Notice  that  for  n  >  0  Eq.(7.20)  reduces  to  a 
two- term  relation  between  i’g(t)  and  Xj^(t)  only. 

Eq.  (7.20)  is  most  easily  solved  in  the  Laplace  domain.  If  we 


introduce 


iwCt-tp) 


(7.22) 


^  (oj) 
n 


I 


dt  e 


then  (7.20)  attains  the  form 

(o)  -  A  +  +  iYn)i((»>)  -  YB(no^i^_ j^(ti))  - 

Here  the  initial  values  '''^(^q)  n  *  0,  1,  2,  ...  appear  as  inhomogeneous 

terms.  The  set  '"j^(tQ)  for  all  n  represents  the  initial  correlations  of  a(t) 

on  t  =  tg,  and  they  connect  the  time-evolution  of  <o(t)>  for  t  >  t^  to  its 
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recent  past.  In  other  words,  Eq.  (7.23)  relates  the  set  ’^^(t)  for  t  >  t^ 

to  the  initial  set  ir  (t«). 

n  0 

Equation  (7.23)  can  be  solved  for  an  arbitrary  initial  set  ^j^(bQ)  by 
24 

standard  techniques,  but  the  solution  is  not  transparent.  In  order  to 
elucidate  the  structure  of  the  solution,  we  assume  a  non-stochastic  initial 
state  o(tQ).  From  Eq.  (7.1)  we  then  find  at  t  »  t^ 

C(p,tQ)  =  P(p)  aCt^)  ,  (7.24) 

and  from  Eq.  (7.14)  we  obtain 


i(p.tQ)  «  o(tQ)  .  (7.25) 

Hence  at  t  *  t^  the  vector  g(p,tQ)  is  independent  of  p,  and  therefore  the 
expansion  coefficients  are  simply 


Then  only  the  recurrence  relation  for  n  *  0  is  inhomogeneous,  and  the 
solution  of  (7.23)  for  *  <5(u))>  is  readily  found  to  be 

<o(u))>  -  - ^ -  o(t-)  .  (7.27) 

u>  -  A  +  iAB  +  R(a)) 


The  effect  of  the  finite  correlation  time,  e.g.,  the  deviation  from  the 
Wiener'L6vy  limit,  is  accounted  for  by  the  operator 


R(u))  ■  yB 


lo" 


0)  -  A  +  iXB  +  liy  +  yB 


2o 


-  yB 

- 2 

u)  -  A  +  iXB^  +  2iy  +  yB  yB 


(7.28) 


which  indeed  vanishes  for  y  -►  0,  X  finite.  In  this  limit,  Eq.  (7.27)  is  the 
Laplace  transform  of  Eq.  (7.8). 

The  explicit  expression  (7.27)  provides  the  exact  solution  for 
i ituations  where  the  initial  state  is  non-stochastic  and  for  cases  where  the 
solution  is  independent  of  the  initial  state.  As  an  example  from  quantum 
optics,  we  mention  that  Eq.  (7.27)  with  o(tQ)  =  1,  A  =  0  and  B  =  1 
represents  the  laser  spectral  profile.  Another  example  pertains  to  the 
long-time  behavior  of  the  solution  <o(t)>.  If  there  is  any  damping  in  the 
system,  which  might  be  caused  by  the  stochastic  f luctutations  itself,  then 
the  solution  for  t  >>  t^  will  become  independent  of  the  initial  state.  If 
we  indicate  by  o  the  solution  <o(t)>  for  t  <®,  then  o  obviously  obeys  the 
equation 

(A  -  iXB^  -  R(0))  5-0  .  (7.29) 


For  the  problem  of  atomic  fluorescence  in  a  strong  laser  field,  this  is  the 
equation  for  the  atomic  steady-state  density  matrix,  which  determines  the 
fluorescence  yield.  There,  the  solution  o  of  Eq.  (7.29)  is  unique. 


VIII.  CONCLUSIONS 

Solving  a  multiplicative  stochastic  process  o(t)  for  its  average  is 
rarely  feasible  by  analytical  methods.  This  is  mainly  due  to  the  finite 
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correlation  time  of  the  driving  process  ^(t),  which  prohibits  the 
factorization  of  the  average  of  a  product  into  the  product  of  the  averages. 
Averages  of  a  functional  of  ^(t)  might  factorize  if  the  process  is  delta- 
correlated.  For  Markov  processes,  however,  we  can  simulate  a  6-correlation 
by  the  introduction  of  the  marginal  average  C(^Q,t)  ■  <6(^(t)-^Q)o-(t)>.  The 
combination  of  the  multiplication  by  dC^Ct)-^^)  and  the  Harkov  property  of 
the  probability  distributions  of  ^(t)  then  gives  rise  to  a  factorization¬ 
like  result  for  the  formal  expression  for  the  average.  Along  the  same  lines 
as  in  a  factorization  assumption,  we  can  now  derive  exact  equations  for 
^(♦Q,t).  In  this  paper  we  have  studied  Eq.  (1.1),  where  we  considered  the 
stochastics  of  ^(t)  to  be  given.  The  equation  of  motion  for  the  marginal 
average  turned  out  to  be  Eq.  (6.4).  The  applicability  of  this  equation  was 
illustrated  by  some  examples. 
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